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Abstract 

This paper addresses three topics concerning the quantization of non-commutative 
field theories (as defined in terms of the Moyal star product involving a constant tensor 
describing the non-commutativity of coordinates in Euclidean space). To start with, 
we discuss the Quantum Action Principle and provide evidence for its validity for non- 
commutative quantum field theories by showing that the equation of motion considered 
as insertion in the generating functional Z'^\j] of connected Green functions makes 
sense. Second, we consider the generalization of the BPHZ renormalization scheme to 
non-commutative field theories and apply it to the case of a self-interacting real sc;alar 
field. Finally, we discuss the renormalizability of various models for a self-interacting 
complex scalar field by using the approach of algebraic renormalization. 
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A Evaluation of expression (50) 



1 Introduction 



It is well known that the laws of standard quantum field theory break down at very short 
distances (at the order of the Planck length) . One of the ideas which have been put forward 
to circumvent this problem is to generalize space-time to non-commutative spaces, the 
simplest one being the flat (Groenewold-)Moyal space - see e.g. [IH2] for an introduction. 
The latter is characterized by the commutation relation 



where ^^'^ = const., i.e. the space-time coordinates are promoted to operators which 
act on a Hilbert space and which fulfill a Heisenberg-type algebra. Unfortunately, quantum 
field theories formulated on such spaces suffer from new types of infrared divergences which 
are related to the ultraviolet divergences of the model and which cannot be removed by the 
introduction of masses. This complication which is referred to as the UV/IR mixing problem 
implies that the familiar models like the </)^-theory are non-renormalizable when written on 
Moyal space [H |5] • In recent years it could be shown that this problem can be overcome 
for the i;^^-theory in Euclidean Moyal space by introducing specific additional terms into 
the action. Thus, the harmonic model introduced by Grosse and Wulkenhaar [6l [7| or the 
1/p^-model devised by Gurau, Magnen, Rivasseau and Tanasa [8] could be proven to be 
renormalizable by using the methods of multi-scale analysiaM 




Several models for gauge field theories on non-commutative space have also been pro- 
posed and discussed [9^16], but so far their renormalizability could not be established. In 
fact, the methods considered for scalar field theories such as the multi-scale analysis cannot 
be applied, or at least not without some serious complications, since these methods break 
the gauge symmetry. For commutative quantum field theories (QFTs), the approach of 
algebraic renormalization \n\ \TE\ is a powerful tool for proving the renormalizability of 
models featuring symmetries. However, the standard formulation of this approach only 
applies to local field theories. Since non-commutative QFTs are non-local, the method of 
algebraic renormalization cannot be used unless it is adapted/generalized to the particular 
non-localities appearing in the non-commutative case. Although this task is far from being 
trivial, the present paper aims to provide some steps towards this direction. 

Our paper is organized as follows. In Section [2| we briefly review the Quantum Ac- 
tion Principle (QAP) in commutative field theory and then we check its validity at the 
lowest order of perturbation theory for a cp^-theovy on Euclidean Moyal space (Section [s]). 
The second topic we address is the application of BPHZ substraction and renormalization 
scheme to non-commutative field theory (Section |4]). The nature of non-localities in non- 
commutative models is analyzed in Section [5] before considering the approach of algebraic 
renormalization to a complex scalar field theory with a rigid U{1) symmetry in Section [6| 

1.1 Brief review of models 

When discussing different aspects of the renormalization of non-commutative QFTs, it is 
useful to have in mind the specific models which have been considered and the results 
which have been established concerning their renormalizability. We therefore conclude the 
introduction with a brief review of models while focussing on real scalar fields on Euclidean 

^In the Grosse- Wulkenhaar case, the initial proof was established by its authors by using the Wilson- 
Polchinski renormalization group approach in a matrix base. A later proof [7j relied on multi-scale analysis. 
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Moyal space in four dimensions. (For other dimensions and for fermionic fields, see for 
instance |19j and references therein.) The case of a complex scalar field will be discussed 
in Section m 

Naive (/)^-theory: The action is given by 



and this model cannot be expected to be renormalizable due to the UV/IR mixing problem. 
Harmonic model for the (/)^-theory: A renormalizable (though not translation invari- 
ant) model for the (/)^-interaction is obtained [HI [7] by adding a harmonic term to the naive 
action ([2]): 

Sgw^ = SoicP] + j d^x Q n\x''<p) * (5^</.)) . (3) 

Here, x^ = 2 {O'^) ^ x'" and > is a dimensionless constant. In the particular case where 
= 1, the action ^ is invariant under the so-called Langmann-Szabo duality. We note 
that the harmonic term admits a geometric interpretation in terms of non- commutative 
scalar curvature [20] . 

Harmonic model for the (/)^-theory: The action reads 

S3[cP]= J d^x(^^d^^d^^+^m^cP^ + ^n\x'^^)*ix^cP) + ^cP*cP*cP^ , (4) 
and is renormalizable [21]. 

(/)^-theory with a 1/p^-term: A translation invariant renormalizable model for the quartic 
self-interaction originates from the inclusion of a non-local counterterm which eliminates 
the most singular part of the IR divergence and thereby overcomes the problematic UV /IR 
mixing problem [8]: 

5'trans.inv. [0] = So[<p] - j d'^X (^(t>^<i^ ■ (5) 

Here, the parameter a is assumed to have the form a = a'/ 9 where a' represents a real 
dimensionless constant and the deformation matrix (0^^^) can be (and is) assumed to have 
the simple block-diagonal form 

/ 1 \ 

-10 
1 
V 0-10/ 

(^^-theory on degenerate Moyal space: Motivated by the work [22], Grosse and Vignes- 
Tourneret [23j studied the </)^-theory on a degenerate Moyal space (i.e. a space for which 
some of the coordinates commute with each other). They found that in this case the 
addition of the harmonic term is not sufficient to establish renormalizability - an extra 
term of the type p (Trc/))^ is required: the action 

Sgvt[4>] =1^^^! \ ct>{x, y) {-/\ + m^ + ^ y^^ ^{x, y) (7) 

~^ ~ff2 j I j ^^-^ '?^(^' y )'?^(^' + ^ j d^X (f) -k (j) -k (f) -k (j) 



with (9 G ]R . (6) 
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is renormalizable. This result is of interest in relation with the attempts to construct 
renormalizable models on non-commutative Minkowski space involving a time coordinate 
which commutes with the spatial coordinates. More generally, terms of the type "product 
of traces" appear to be quite natural for theories on non-commutative spaces |24j . 

2 QAP for commutative QFTs 

The Quantum Action Principle (QAP), i.e. the renormalized version of Schwinger's action 
principle, plays an important role in the study of field theories which are invariant under 
transformations that depend non-linearly on the fields, e.g. non-Abelian gauge field the- 
ories. These symmetries are expressed by Ward identities which must be satisfied by the 
Green functions of the theory. The situation is well understood for QFTs on commutative 
space and details concerning power counting renormalizable QFTs may for example be 
found in the references [TTl [THl . Quite generally the QAP describes the result of the 
insertion of a composite operator into Green functions in the renormalized theory. There 
are three versions of the QAP corresponding to the different types of insertions to be con- 
sidered: equations of motion (e.o.m.), local field polynomials and derivatives with respect 
to parameters. In the following we briefly recall these three versions by considering the 
example of the i;^i4-theory, i.e. the power counting renormalizable scalar field theory on 
commutative Minkowski space-time defined by the classical action 

= So[<P] + SinM ^ J d^x i-^cPlC^c^) + J d^x (-^ 0^) , (8) 

where 

}C^ = n^ + m\ (9) 

The fundamental quantities in quantum theory are the Green functions, in particular the 
connected Green functions: the latter are collected in a generating functional Z'^[j] = 
—ihln Z[j] from which they can be extracted by differentiating with respect to the external 
source j, 

where the (j){xi) are Heisenberg field operators and h = 1, see for instance reference [26]. 
Insertion of e.o.m. 

Starting from the generating functional Z[j] = exp{iZ^[j]) of Green functions, a formal 
calculation leads to an equation representing the e.o.m. of the theory expressed in functional 
form: 



^^Z<^[j]+j{x) = 0. (11) 

i Sj 



64>{x 

By using any of the known renormalization schemes (BPHZ, dimensional or analytic regu- 



larization, . . . ), one can prove that an equation of the form (11) holds in the renormalized 



theory, i.e. the statement which is referred to as the first version of the QAP: 

O{x)-Z'[j]+j{x) = 0. (12) 
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Here, 0{x) ■ Z^[j] denotes the insertion of the local composite operator 0{x) into the 
renormalized connected Green functions, i.e. 

(o|r(^(xi) . . . 4>{xn)\o)c {o\Tdix)4>ixi) . . . H^nMc , (13) 

and at the lowest order in the ^-expansion 0{x) is the local classical field polynomial 
6S/5(j){x) (whose vanishing is tantamount to the classical e.o.m.): 

d{x) = ^ + Oih) = -fCjix) - ^ 4>\x) + Oih) . (14) 

Furthermore, the dimension of 0{x) is bounded by 4 — dim (p. 



For the free theory (i.e. for A = 0), Eqn. (12) is nothing but the free e.o.m. in the 
presence of the source j. 

Insertion of field polynomials 

The second version of the QAP which can be proven, consists of the insertion of Q{x) 



' 5<t>(x) 

0{h) where Q{x) is a classical field polynomial. 



Insertion of derivatives with respect to parameters 

as 



The third version comprises of the insertion of §y + 0{h) where A denotes a parameter 



appearing in S. 



Some consequences 



The first version of the QAP as given by equation (12) states that the functional F[j] 
appearing on the Ihs of this equation identically vanishes or equivalently, that all of its 
Taylor coefficients vanish, i.e. 



6^F\j] 



6j{xi) . ..5j{Xn) 



= 0, for n = 0,1, 

j=0 



Each of the latter identities can be viewed as a perturbative expansion in the coupling 
constant A. In the following, we will explicitly derive the relations that one obtains for 
n = 1 at the lowest orders in A. These identities are ensured to hold as mere consequences 
of the QAP whose validity is well established. In Section[3j we will then verify explicitly the 
validity of the analogous consequences for non-commutative (/)^-theory at next-to-leading 
order in perturbation theory, thereby providing evidence for the validity of the QAP in the 
non-commutative setting. 

By differentiating relation |l2| with respect to j{y) at j = and using = S^'^'^{x—y), 
one concludes that for the free scalar field theory 

ij(^)(:E -y) = (O|rd(x)0(2/)|O)e = -/C,.(O|r0(x)(A(y)|O) , (15) 

i.e. the defining equation for the propagator whose solution reads 

(O|r^(x)0(2/)|O) = iG(x -y)= hm / ^e-'^'-^-v) . (16) 

e^o+ J (27rj^ /c^ — -|- le 
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For the interacting theory, differentiation of relation (12) with respect to at j = 
yields 

A 



|0), + i5(4)(a;-y) =0. 



We may then apply the Gell-Mann-Low formulsj^ 

(OITexp (i Jd'^zCint) |0)o 

and make use of Wick contractions for free asymptotic fields: 

I 1 

r(0(x)0(y)) = : Hx)Hy) ■■ +Hx)ky) , 
I — I 

with Hx)ky) = {0\T{4>{x)4>iy)mo 1 = iG{x -y)l. 
To first order in the coupling constant, expression ( |17[ ) then reads 

iA /" ,^ ':a , ^\ , s A 



(17) 



(18) 



(19) 



JCMx) 1 



\0)o^, + i5^^)(x-y) = 0, (20) 



where the subscript 0, c means that we only consider the connected parts (the non-connected 



contributions being compensated by the denominator in (18)). By using Wick contractions 



expression (20) becomes 



A 



I 1 



I 1 I 1 



ICx(l)ix)(j){y) - ilCx- / d z(f){x)(l){z)(j){z)(j){z)(j){z)(j){y) 



+ 



A 



{x)kx)kx)(p{y) + iS'-^^x -y) = 0. 



(21) 



At order in A, we recover relation ( |15[ ). At first order in A, we obtain the result 

/C^. j ct^z G{x - z)G{0)G{z -y) + G{0)G{x - y) = , (22) 
which may be interpreted graphically by: 





At order A^, the expansion (17) yields 
(0|T 



)Cj{x)h^f id^zi id^Z2^\zi)^\zMy) 



2' 4! 



+ ^0'(^)(^) ld'zk{z)cp{y) 



|0)o,c = 0. 



(23) 



This relation as well as those obtained at higher order in A can be explicitly worked out 
along the lines indicated above. 



■^Note, that the Ihs involves Heisenberg field operators, whereas the rhs concerns asymptotic free fields 
as indicated by the subscript 0. 
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Non-linear field variations. Now let us consider field variations which we denote by 



p{4>{x)) =: P{x). Eqn. (12) is then replaced by p^J 



(24) 



where p denotes an external source which is coupled linearly to p{4>). Varying twice with 
respect to j and setting the sources to zero afterwards, leads to the following relation for 
graphs with two external legs [271 EH] • 



(0|r d^xP{x) 



5S[ 



6(f>{x 



-<Pixi)<Pix2)\0)c + {0\tP{xi)(I){x2)\0)c + (O|rP(x2)0(xi)|O), 



0. 



(25) 



Euclidean case. Before closing this section, we should mention that the QAP is also 
valid for field theories on Euclidean space (see e.g. reference ^7\), the obvious difference 
with the Minkowskian case being factors of iti. More specifically, in Euclidean space the 
generating functional for the connected Green functions is given by 



Z'[j] = -In Z[j], 

and the Gell-Mann-Low formula may be expressed as 



Z\j]=Me 



-Sin 



(26) 



(27) 



3 QAP for non- commutative scalar quantum field theory 

We will now "switch on" the non-commutativity, i.e. consider the non-vanishing commu- 
tator of Eqn. ([T]). In order to simplify the calculations, we employ the Weyl quantization 
map W which enables us to implement the non-commutativity of space in terms of a de- 
formed product for the functions defined on this space, namely the (Groenewold-)Moyal 
star product (see e.g. reference m for an introduction): 



Since 



ei^-929^;/(x)5(y) 



d^x{fi-k----kfn){x), 



x=y 



(28) 



(29) 



the integral of a star product of functions shares the property of the trace of being invariant 
under cyclic permutations. 

For the discussion of the naive (^4-theory on non-commutative Minkowski space, we 
will follow the conventions of Micu and Sheikh- Jabbari [29j, and consequently we choose 
gO* = in order to avoid difficulties with the time ordering [SUj- The action for the naive 
(/)^-theory reads 



S[cP,X] 



(fx 



{d^(l)*df,(l)-m^(t)' 



A 
4! 



< -k -k -k t 



(30) 
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and the property (29) of the star product allows us always to drop one star. The Feynman 



rules for this model are given by 



Jj2 _ _j_ 




^4 



iV{ki,k2,k3,k4) 

^(2^)V^) iki+k2 + k3 + ki 



cos ( ^ A;iA;2 j cos (^k^k^ 



cos ( ^ kik^ij cos f ^ k2k4j + ^1^4 j cos f ^ ^2^:3 



where we have introduced the notation 



(31) 



(32) 



In the following considerations concerning QFT, we will omit the hats on the Heisenberg 
field operators in order to simplify the notation. 



3.1 Verifying the QAP for the e.o.m. 

As indicated in the previous section, we would like to provide evidence for the validity of 
the QAP in non-commutative field theory. In order to do so, we assume that the latter 



has the standard form, i.e. F[j] = where F[j] is given by the Ihs of equation (12), hence 

5F 



Sj{x) 



= 0. In the sequel, we explicitly verify this relation to lowest order in perturbation 

j=0 

theory for the nai've (/)^-theory in Minkowski space. A complete proof of the QAP (which is 
beyond the scope of the present work) requires both the consideration of all orders of the 
Taylor series expansion and a focus on renormalizable models in Euclidean space. 



The e.o.m. associated to the action (30) is determined by 



0{x)^^-^^ = -)C^^{x)-^^*^{x), with <A^3^ (33) 
0(p[x) 3! 

The QAP for non-commutative quantum field theory is formulated as 



O{x)-Z'^[j]+j{x) = 0. (34) 
Differentiation with respect to j{y) at j = yields 

|0)c + i5(')(x-y) = 0. (35) 

In order to evaluate this expression perturbatively using the Gell-Mann-Low formuls[^one 
needs to consider the Wick contractions in the non-commutative setting. For convenience, 

^At this point we assume that the GeU-Mann-Low formula holds in the non-commutative setting. The 
computation carried out in the remainder of this section indirectly verifies this assumption up to order A. 



lC,cP{x)<Piy) + -r\x)<P{y) 
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we use their momentum space representation as following from Eqn. (19) and Eqn. (16): 



I 1 

{x)<j){y) = lim 



-ik{x-y) 



1 



0+y (27r)^ 



This implies the identification 



I 1 



(27r)4y (27r)4 
(2vr)^i 



^2 _ ^2 _|_ 

I 1 



-ikix—\k2y 



kf — m? + ie 



(36) 



(37) 



Let us now return to Eqn. (35). As in the commutative case (see Eqn. (20)), we obtain up 

A 



to first order in the coupling constant A 

iA 
4! 



(o|r 



lCMx){l 



d:'z<P^\z)](t^{y) + -^*'\x)<l,{y) 



|0)o,c + i^^'Hx - y) = . (38) 



At the order zero in A, this relation represents the defining relation for the propagator. Let 
us scrutinize the terms in Eqn. (38) which are of order one in A: the first of these terms 
reads as follows in momentum space 



4! 'J {2ttY J (27r)4 



k\h 



cos 



IK2 



cos 




+ COS 



kik4 



COS 



^2^3 



(2^)4 5(^)^fc, |0) 



0,c • 



(39) 



.1=1 



We can now employ Wick contractions in momentum space and obtain 



-iA 



/C, 



(270^(2^ 




(40) 



15(4) (gi + ki) 15(4) (^2 + ks) 15(4) (A;4 + (?2) 



qf^ — rn? + ie fcg — m? + ie ^2 ~ ^ + 



k-\ h 



cos 



1«'2 



COS 



+ COS 



k^ k. 



1^3 



k2k4 



COS 



+ COS 



k-i k, 



IK4 



COS 



k2k3 



and, after some further simplifications, 



-iA 

~6~ 



d% d% 5(4)(gi+g2) 



(2^)4 (27r)4 g2 _ ^2 + 



X 



COS 



qik2 



COS 



^^2^2 



+ COS 



A f d% 



2 

iqi(x-y) 



(271)^1 gf - m2 + ie 



cos 



d'^ko 



^292 



A;^ — + ie 



+ cos 



qiq2 



(27r)4 A:2 - m2 + ie 



cos(giA;2) + 2 



(41) 



8 



This is the well-known tadpole contribution with planar and non-planar parts, respectively. 



To conclude, we consider the second term of order A in Eqn. (38). In momentum space, 
we have 



hence we obtain 



(2^) 



12 



A 
3! 



{0\T<p*'ix)<Piy)\0)o,, 



A f d^q 



I 



(27r)4 A;2 - m2 ie 



6 J (27r)4 -m'^ + ie 



cos{qk) + 2 



- 5 fci A;2 - i fci fcs - i fc2 fcs 



-iq{x-y) 



(42) 



(43) 



This expression cancels (41), and thus we have demonstrated the validity of the QAP for 
the e.o.m. to first order. 

3.2 Verifying the QAP for a non-linear field variation 



Here we study the non-commutative analog of Eqn. (25) up to one-loop level for the (j)*^- 
theory by considering the example 



(44) 



Using the Gell-Mann-Low formula (|18|) and the lines of arguments presented in the previous 

0(xi)0(x2)|O)c 



subsection (cf. Eqns. (25), (33)), one arrives at the relation 

iA /" ,^ , . ^ \ A 
4! 

4! 



= (0|r y d^x(/)*3(x)* 
-{0\T^*\xi)cl){x2)(l 



d''zct>^\z))\0). 



(o|r(/>^=^(x2)0(xi)(i 



iA 
4! 



d^Zc/>'^^(z) |0)e. 



(45) 



In order to check explicitly that this relation actually holds to order A, all possible Wick 
contractions leading to connected graphs have to be considered. When working in x-space, 
particular care must be exercised concerning the order of the resulting propagators due to 
the star products, for instance 

I 1 

(j){x) ★ (j){x) * (f){x) 4>{y) = 4>{x) -k iG{x — y)-k (j){x) . (46) 



Keeping this fact in mind, the explicit check of relation Eqn. ( |45[ ) is straightforward, albeit 
tedious. 



4 BPHZ for non-commutative field theories 

We recall \n\ [26] that the BPHZ method is a recursive substraction scheme which does 
not refer to or use any particular regularization procedure: the renormalized Feynman 
graphs are obtained by subtracting out the divergent parts of the integrands for all graphs 
describing divergent integrals. These subtractions are equivalent to the addition of local, 
regularization dependent counterterms to the Lagrangian, but there is no need to introduce 
these counterterms explicitly. 
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4.1 General considerations 



In non-commutative QFTs we encounter the problem that the usual BPHZ momentum 
space subtractions do not eliminate the infrared divergences which are due to UV/IR mix- 
ing. This may be illustrated by the example of the non-commutative (/)|-theory in Euclidean 
space whose Feynman rules are given by 

G{p) = — 2 ' 



^(Pl,P2,P3,P4) = -Xf{pi,9) 

_ -A 



COS (2f^) cos (ef^) + COS (ef^) cos (2^) + COS (ef^) COS (2f^) 



(47) 



For the non-planar part of the one-loop four-point lP/-graph, the BPHZ-subtraction in 
momentum space yields an integral of the following type: 



4 I cos{kp) 1 



"^^^^ ~ J ''\{{k+p)'^+m2){k^ + m^) (k^ + m^y)' ^^^^ 
Here, p denotes the total external momentum, p^ = O^^p'^, and the second term in the 



integral (48) represents the subtraction of the integrand evaluated at p = 0. (For later 
reference we note that pp = p^O^j^^p'^ = 0.) In commutative (/)|-theory, the phase factor 
cos(A;p) is absent and the second term then cancels the logarithmic UV singularity of the 
first integral leaving just finite terms. In non-commutative (/>|-theory, the presence of the 
phase factor ensures that the first term is UV finite, but it now suffers from an infrared 
singularity in the external momentum which is not canceled by the second term - worse 
still, the second term introduces an additional UV singularity. 

As a second example, consider the BPHZ-subtraction for the combined contributions 
of planar and non-planar parts: 

7 ( \ - /"^4l ( ^ + ^ cos{kp) A + B \ 

Jtotip) = jd k \^((^+p)2+^2)(fc2 + ^2) - (A;2 + ,„2)2 J ' (^9) 

Clearly, the second term (term proportional to A + B) overcompensates the UV divergence 
of the very first integral (term proportional to A) while the infrared divergence remains 
present. 

From these examples we conclude that the usual BPHZ-strategy of subtracting terms 
with p = from the integrand does not work for non-commutative theories. Here, we 
suggest the following procedure for the BPHZ momentum space subtraction in non- commu- 
tative theories: consider p and p as independent variables (though satisfying pp = 0) when 
applying the subtraction rules of the BPHZ scheme, i.e. subtract from the integrand terms 



with p = while maintaining the phase factors: instead of the subtraction (48) we thus 
consider 

J(v)= !d^k( "^"'^^^^ cos(fe^\ 

]'"'\{{k+pf+m?){k^ + m'^) (^+^2)2^- ^^"^ 

By proceeding in this way, the logarithmic IR divergence of the first term is canceled by 
the second. 
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Expression ( 50 ) is explicitly evaluated in Appendix [A] leading to an integral over a 
modified Bessel function Ki{z) of second kind: 



Jip) = ttV 



1 (1 - 2e)(e - 1) i^i vi^^ - - 
di , ^ 



(51) 



This expression is perfectly finite since the series expansion Ki[z) = - + i^zlnz + 0{z) 
implies that we have for small values of p: 

y/p2 _|_ 4j^2 _|_ y/p2 



J{p) 



-IT 



Am? 



In 



y/p2 _|_ 4j^2 



+ 27r2 + 0(p2 



(52) 



This expression is regular in the limits — t- and/or m — t- 0. 

Next, we consider the planar part of the momentum space integral which is well known 
in the literature since it also appears in the commutative (/>4-theory [IT]: 



Jpi{p) 



(53) 



{{k + p)"^ + rin?){k'^ + rn^) [k"^ + m?)"^ ^ 

The application of the procedure described above does not require any extra calculation: 
By setting ry = in equation (103), we end up with 



/"l /"OO 

/ di{l-2i){i-l) / dAe-M™'-(€-i)?p'] . 



(54) 



After carrying out the integrations, we recover the well-known expression for the planar 
integral: 



Jpi{p) 



-TT 



1 H ?r- In 



p. 



y/p2 4m2 + a/p^ 
\/p'^ + 4m^ — Y^P^ 



(55) 



Comparing this result with expression (52) which holds for small values of p, we conclude 
that our procedure for applying the BPHZ subtraction leads to differences between the 
planar and non-planar parts which only appear in the higher-order terms 0{p^). 

A general one-loop integrand I{pi,pi,k) with superficial degree of divergence n (upon 
integration over k) should be replaced by 



f,f{Pi,Pi 



R{pi,Pi, k) = {l- tp) I{pi,Pi, k) 
{t;fKPi,Pr) ■.= f{0,pi) + Y.Pji 

j V 

_|_ jL \ " A»l ( ^ 

^ n\ ^ ^'il ■■■Pin U Ml 



d_ 
dp' 



P^=0 

d 



+ 



dp^, 

Jn 



:i^f{Pi,Pi) 



Pi=0 



(56) 



Note that only the external momenta pi are involved in this subtraction scheme and that 
Pi ^ at all stages so as to ensure convergence of integrals as in the commutative case. 
This also means, that one must exercise care in rewriting phase factors as kOpi = kpi rather 
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than —kpi. Also note that phases not involving k will not regularize the integrals and hence 
not lead to UV/IR mixing. There is, however, some ambiguity concerning our treatment 
of phases involving piOpj. We hence suggest to rewrite these as pidpj = \{piPj — PiPj)- 

The strategy outlined above, which amounts to a simultaneous subtraction of UV and 
IR divergences, allows us to overcome the UV/IR mixing problem at one loop order. In 
the next subsection, we will show for the example of non-commutative 0|-theory that this 
procedure amounts to a one-loop order redefinition of the parameters appearing in the 



classical action. In Subsection 4.3, we will comment on higher-order loops. 



4.2 One-loop renormalization of non-commutative 04-theory 

We now apply the subtraction scheme presented in the previous subsection to the non- 
commutative (j^l-theory. Thus, we consider the following generalized action in Euclidean 
space: 



S\ 



1 



(1 + ^)5^ 



(57) 



Here, the parameter A is dimensionless, B has the same dimension as and C the same 
dimension as A. The Fourier transforms of the parameters A, B and C, which we denote 
by a, h and c, potentially depend on the momentum and on 0. The dependence of c on the 
momentum is motivated by the fact that the star product leads to a momentum dependent 
coupling. As usual in the BPHZ-scheme, m and A are already the one-loop renormalized 
mass and coupling, respectively. We consider a, h and c to be of order h (i.e. one-loop 
expressions), and for h = 1 the Feynman rules in momentum space are given by Eqn. (47). 
In the following calculation it will be essential to keep p^ ^ d'^P^, since we treat p and 
p as independent variables in the BPHZ subtraction scheme as discussed in the previous 
subsection. 

Including its one-loop correction, the vertex T^^^{pi,p2,P3,P4:) = T^'^^{pi) is given by 
(see Eqn. (|47]) and references (29l [3T] ) 



r(^)(p,) = (A + c(p„e))/(K,0) 

A2 f d^k f 1 

+ 



9 J {2it)^ \ ik+pi+p2y + m? k"^ + m? 



1 + 



gifcPi _|_ gi's(pi+P2) _j_ 1 ^ gifc(pi+Pi) 
i=l i=3,4 



(58) 



Thus, the vertex correction consists of integrals of the same type as those discussed in the 
previous subsection. 



We now add and subtract the same term to expression ( 58 ) : 

T(*\p,) = {x + c{pi,e))f{pi,e) 




1 



1 



{k + pi + p2)^ + m"^ k^ + m"^ {k'^ + 

4 



' ^ i=l 



1 



9 J (2^)4 (A:2 + m2)2 



1 + 



+ 



i=l 



(59) 
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The second line of this expression is finite - let us call it A'{pi,6). For small external 
momenta pi, one has f{pi, ^) ~ 1 and we get [29] 



1 



(2^)4 (fc2+m2)2 



6(47r) 



y^ln(?T7-^p^^) + ^ln(m^(pi + Pi)"^) 



.i=l 



(60) 



i.e. we have logarithmic IR divergences. 

Following the lines of the BPHZ scheme, we now split c according to 

c{pi, 0) = C + c'ipi, 0) + Coo, 

where c' is an arbitrary finite constant and 



(61) 



6(47r 



i=l 

2A2 f d^k 



^ln(m^]3^) +^ln(m^(pi + Pi)"^) 
1 



1=2 



9 y (2^)4(yfc2 + m2)2 • 



(62) 



The function F^^^ resulting from the decomposition (61) is denoted by ri'^^: for small pi it 
thus reads 

F(4)(p,) = A + c' + A'(pi,e), (63) 
where the constant c' is to be fixed at a convenient symmetry point (s.p.) by the condition 



ri'Hp^ 



is.p. 



A. 



I.e. c 



-A'l 



is.p. 



(64) 



This defines the physical coupling at the symmetry point and at one- loop order. 

We now turn to the case of generic (i.e. not necessarily small) external momenta. From 
( 59 ) we then conclude that 



r(^\pi) = {x + c{p.,e))f{p,,e) + A' + 



2A2 f d^k 



a2 r d^k 



9 J (27r)4 (A;2 + m?) 
A2 f d^k 



9/7 (2^)4(A;2 + m2)2Ll2Z^« ^■■■) 



+ ... 



+ ... 



9/ J (2^)4 {k'^+m?) 



(65) 



henceforth 



F(4)(p,) = (A + c')/(p.,e) + A'(p, 



Pi 



Is.p. 



s.p. 



(66) 



The symmetry point is defined in terms of the Mandelstam variables s,t,u and some 
normalization mass /i: 

s = {pi+P2f, t = {pi+P^f, u = {pi+p4f, 



2 2 w 
Pi = \fl.. 



PiPj 



(67) 
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Next, we consider the two-point function at the one- loop level, r(2)(p,-p) = r(2)(p2)^ 
which is given by 

r(^)(p^) = (1 + a)p- + ™^ + Hp,») - I l^'-i^ ■ 

Since the latter integral only depends on p and not on p, we consider a = and we split b 
according to 

b{p,e) = b' + b"{p,9) + boo, 

A f df^k cos{kp) A 



6 J (27r)4 /fc2 + ^2 247r2 W p2 



^ + ^ [In (m2pV4) + 27i, - 1] + O {p'] , (69) 



6o 



247r2p2 967r2 
A /• d^fc 2 



6 7 (27r)4A;2 + m2 ' 

where 6' is a finite constant. The one-loop renormalized two-point function hence reads 

r(2)(p2) =p2^^2^^/^ 

supplemented by the normalization conditions 

^PiP%.=.^2=0, i.e. b' = 0, (71) 



2 = 1. 

The calculation and reasoning presented in this subsection should also apply to gauge 
theories (in particular to the translation invariant models introduced in reference |13j). 
although the following complications are to be expected in this case. 

In Minkowski space, as long as time commutes with the spatial coordinates (i.e. 0^^ = 0), 
the on-shell condition p2 = q does not mean p^ = ^ since there is no time-component in 
p. Hence, IR issues only arise in integrals, but not in connection with the normalization 
conditions which are to be considered in the actual renormalization process. 

In Euclidean space, the on-shell condition p^ = Q may imply p2 = q depending on the 
structure and rank of {9^^). Hence it might be preferable to consider Min(p2 -|- 1/p^) as 
the renormalization point instead of = 0. This would apply also to the normalization 
conditions fixing the physical parameters. 

Before proceeding further, it is worthwhile to come back once more to the well-known 
non-planar contribution to the propagator as described by the integral 



^//^^ 0) - ^ /" cos(A:^) 



6 J (27r)4 fc2 + ^: 



2 



see expression ( 69 ) . For small values of fP' , this integral exhibits a quadratic and a logarith- 
mic IR divergence at fP' = 0. These IR divergences, which are tied to the UV divergences, 
represent the most basic manifestation of the UV/IR mixing problem appearing in non- 
commutative QFT. Multiple insertions of this non-planar diagram into a non-planar tadpole 
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graph lead to IR problems for higher order graphs [H [291 132] • These problems motivated 
the authors of reference [8] to introduce the most singular (i.e. quadratic) divergence into 
the classical Lagrangian by means of the non-local term 

J (fx (0^0) , with □ = 9^9^ = , (72) 

where a' represents a real dimensionless constant. The renormalization of the resulting 
model was established using multi-scale analysis in reference [5] and the one-loop renor- 
malization was discussed more specifically in reference [32]. The crucial point is that the 
modified propagator G = {p^ + m? + a'"^ /fP')~^ has a damping behavior for vanishing mo- 
mentum (i.e. limp^o G{p) = 0) which ultimately allows to cure the IR problems at any 
order. 

In the BPHZ approach considered in this section, the inclusion of the non-local term 



(72) into the classical action modifies the quadratic part of the functional (57) according 



to 



^ (1 + A)d^'<^d^^ + \{m'' + B)<t? ^ ^ (1 + A)d''cPd^cl) + I 0(m2 + £ + 5)<^ , 
^ ^ 2i 2i I I 



the normalization condition ( 71 ) becoming 



Hence, the one-loop renormalized two-point function (68) is modified according to 

where b(p, 6) can now be interpreted as the renormalization of the new parameter a' — see 
reference |32j . 

4.3 Beyond one loop 

Our momentum space subtraction scheme allowed us to discard all divergences at one-loop 
order. Potential complications arise in higher-loop graphs: 

• For all divergent subgraphs of a given graph, we need an unambiguous distinction 
between internal and external momenta since the latter are treated as passive variables 
in our subtraction scheme. A point which is related to this issue is the need to 
ensure that Zimmermann's forest formula still allows for a consistent treatment of all 
divergences. 

• In order to establish renormalizability, it has to be shown that the ambiguities in 
the subtraction procedure result in a finite number of counterterms (redefinition of a 
finite number of parameters determining the Lagrangian) . 

Further investigation of these problems and in particular the tackling of two-loop graphs 
is currently under study [3^. Here, we only emphasize the following point. Since the 



'naive" 04-theory described at the one-loop level by the action (57) is not renormalizable 
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it appears to be clear that the presented subtraction of divergences leads at some stage 
to the inclusion of some extra term in the classical Lagrangian. In the next section, we 
will discuss the nature of non-local terms which can appear in the framework of non- 
commutative QFT. These arguments lead to the conclusion that the only non-local term 
which can be included in a translation invariant model for a real scalar field with quartic 
self- interaction is the one given by expression (72). 



5 Nature of non-localities in non- commutative models 

The approach of algebraic renormalization has been developed for local field theories. One 
would expect that due to the loss of locality we have an infinite number of possible coun- 
terterms for any such theory. However, this is not really the case since the non-locality 
only originates from the star product and only leads to IR divergent non-local countert- 
erms in the effective action. This is a consequence of the UV/IR mixing and thereby these 
problematic counterterms can only have an IR divergence of degree equal or less than the 
superficial (maximal) degree of UV divergence d^. Furthermore, these terms are generated 
by phases depending on the combination p^^ = O^^pu where p is some external momentum. 
Hence, IR divergent non-local terms N{p) behave like 

lim Nip) ~ -4^ , with < s < ^ . (74) 



5.1 Translation invariant models 

Let us consider (/)4-theory in commutative space as an example for a power counting renor- 
malizable theory. Its superficial degree of UV divergence is given hj d-y = 4 — E where E 
denotes the number of external legs. If we generalize this model to non-commutative space 
by replacing point-wise products of fields by star products, non-local counterterms arise. 



But the degree of IR divergence, as introduced in Eqn. (74), can be at most s = 2 — E/2 
due to power counting. Moreover, if we use Schwinger's exponential parametrization for 
the integrals (see e.g. Eqns. (102), ( 103] )), we can use the phase factors involving p to com- 



plete the squares in Gaussian-type integrals: for small values of p, the resulting expressions 
behave as 1/p^ where s is a multiple of D/A [D being the space-time dimension). Thus, a 
combination (p^)'-'^"^/^-* x (counterterm) should be local. 

This argument motivates us to generalize the algebraic renormalization procedure by 
allowing non-local counterterms which become local upon multiplication by {p^Y with 
s = 2 — ^ while excluding non-local terms of a different nature. 

Example 1 Consider the real scalar field theory action in Euclidean space 

S = jd'x Q {d,ct>){d^<l^) + \ m^'^ + A<^*4^ _ (75) 

It is translation invariant as well as invariant under the transformation (p ^ —cf). In com- 
mutative space, i.e. for the ordinary point-wise product, the action already includes all 
terms which are allowed for a local power- counting renormalizable theory in four dimen- 
sions. However, in non-commutative space one more additional term is allowed according 



to our considerations above, namely (72). This term is invariant under translations as well 



as under the transformation (j) — )• —(p. Furthermore, it has mass dimension 4 (which is 
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the maximal value allowed for the 4- dimensional model under consideration), and it has an 
IR divergence of order which is the maximum degree allowed by power counting since 
s = 1 for a counterterm with two fields (p. 

Example 2 Now consider the complex scalar field action in Euclidean spac^ 

/( o'^ A2 \ 

d^X id^ipd^if + 'f'^y-' + ^'^'P^ + — (p-k(p-kip-kip+— (f-kip-kif-kipj . (76) 

This action is again translation invariant, invariant under if — )• —ip, (p — )• —ip, invariant 
under the charge conjugation transformation C : ip ip. Furthermore, there exist a 
continuous internal symmetry transformation leaving the action ( 76 ) invariant, namely the 
global U{1) transformation 

6s(p = iecp , de^p = -\e(p , (77) 

the latter leading to the classical Ward identity 

Wr(°) = , with W = jd^x (^p^ - . (78) 

Note that the only possible non-local term which is allowed by power counting and compatible 
with the symmetries, namely (p^'p has already been included in the action (76). 

Assuming the quantum action principle to hold in the non- commutative case and in- 
cluding non-local terms of the nature described above, we can follow the standard procedure 
for studying the breaking of symmetry by radiative corrections (see reference {llj and Eqns. 
(92)- (95) below). We then find the following basis of integrated insertions of dimension 
less or equal to four: 

yy;r("-l) = ^i^A + 0(^+1) , 

A : J d'^x [p^ — ip'^) , J d^x {d^pd^p — d^pd^p) , j d'^x (^V^f ~ ' 

'd^x {p*^ - p*^) , j d^x {p*'^p - p'^'^p) . (79) 
Thus, the basis of admissible terms includes one non-local expression. 



Remark. The additional non-local terms which are bilinear in the fields will lead to 
damping properties of the field propagators in the IR regime. These terms should take care 
of the IR problems in higher loop integrals (see next section). However, it must be shown 
explicitly that they are sufficient to render the model under consideration renormalizable. 



5.2 Non-translation invariant models 

As remarked in subsection |1.1[ non-local terms of the type "product of traces", e.g. 
^(Tr(/))^ or (Trc/)^)^, appear to be natural on non-commutative spaces [23J. In particular, 
a term of this type ensures the renormalizability of the harmonic model for the (/)^-theory 
on degenerate Euclidean Moyal space - see Eqn. ([T]). Yet, it would be preferable to have a 
more precise characterization of the generic non-local terms which are admissible on Moyal 
space. 

*Note that only the vertex proportional to A2 exhibits UV/IR mixing - at least at the one-loop level - 
as was shown in references |34j . 
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6 Algebraic renormalization of NCQFTs 



The quantization of Lagrangian models with continuous (rigid or local) symmetries can con- 
veniently be addressed using the approach of algebraic renormalization fT7] . In particular, 
the BPHZ subtraction scheme can be applied to these models. 



6.1 The self- interacting complex scalar field with a rigid f/(l) symmetry 

In this section we consider models in four dimensional non-commutative Euclidean space 
for a complex scalar field 99 with a self-interaction and we try to apply the methods of 
algebraic renormalization. More precisely, we consider models with a rigid U{1) symmetry, 
i.e. invariance under the infinitesimal phase transformations 



6e(p = iep , 



6eP 



-leif, 



(80) 



where e denotes a constant infinitesimal real parameter. 

As noted some time ago by Aref'eva et al. |34j . the star product allows for two different 
C/(l)-invariant quartic interactions: 



(fx 



Al 



A2 



{d^p) {d^Lf) + m (pip+ — (p-kipip-kip + — if-kipif-kip 



(81) 



Here, Ai and A2 are two independent coupling constants and we have used the properties 



of the star product to drop one star in the quartic terms. The action (81), which will 



be referred to as the "naive" model, is invariant under translations, under the rigid U{1) 



transformations (80), as well as under the discrete transformations of charge conjugation 
C : (p ■v^ (p and of reflection ip — )• —ip, (p — )• —ip. The corresponding Feynman rules have 
been given in references pi]. Since the bilinear terms of the action do not involve the star 
product, the propagator in momentum space is the same as in the commutative theory, i.e. 
{jp' + m?)~^. In momentum space the interaction terms read a^ 



/ 



d'^x V{p, p) 



Al cos ^— + — J + A2 COS ^— J COS ^— J 



(82) 



6.2 Self-interaction T<r (y9 



To start with, we consider the naive model (81) with A2 = 0, i.e. for a quartic term with 
"alternating" fields: 



Sal[p,p\ = I d'^X 



(d^p) (d^p) + m?pp + ^(p-kp-kpkp 



(83) 



It was shown in references |34j that this model is renormalizable at the one- loop level. 
Indeed the considered interaction does not lead to a UV/IR mixing - at least at the one- 



loop level - since the interaction (82) with A2 = does not yield an infrared problematic 



^Note that our convention for the star product differs by a factor 1/2 from the one used in ref. 
hence the factors 1/2 in the cosines below. 
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tadpole contribution at the one-loop level. Later it was argued [35J that the model is 
renormalizable to all orders. 



The following extension of the action (83) has also been studied in the literature: 



Sgwlsz if, f] 



{D^Lp) -k {D^ip) + m?LpLp + fi^ {x^Lp) -k (x^ip) + ^ (f -k if i< ip i< (f 



(84) 

Here, D^ip = d^ip—iax^-kip can be viewed as the covariant derivative describing the coupling 
to a constant "magnetic field" and the term in Q"^ is the complex version of the Grosse- 
Wulkenhaar harmonic term. The model (84) has been proven to be renormalizable to all 
orders |36] and the model with i7 = and q 7^ is known as the Langmann-Szabo-Zarembo 
(LSZ) model, this model being exactly solvable for a = 1 |37] . 



6.3 Extended action for the complex scalar field 



As discussed in references |34j, the naive model (81) is not renormalizable for arbitrary 
values of Ai and A2, although it is in the particular cases A2 = and Ai = A2. In order to 
have a renormalizable interaction for a generic value A2 7^ 0, we include an extra non-local 
term into the action of the type considered earliei|^ for a real scalar field [3 [T3t [32] : 

Ai _ _ A2 



r(o)[(^,^] 



{df,ip)(d^ip) + ip — ip + n?ipip + ^iP'kipi.if'kip+^ 



(85) 



Here, the parameter a is assumed to have the form a = a' /6 where a' represents a real 
dimensionless constant. The propagator in momentum space reads 

1 



G(A:) = - — ^' (86) 

+rin? + '^ 

and its "damping" behavior [8] for vanishing momentum (i.e. lim G{k) = 0) allows us to 

avoid potential IR divergences in higher loop graphs |32j . The vertex in momentum space 
can be read off from Eqn. (82). 

The one-loop corrections to the propagator (including a symmetry factor of 1/2) are 
given by 

n(rt = ^ / (A.+ A.)+ A.cos(>.p) ^ ^ _ 

4 J (27r)'^ 1^2 j_^2 

where HP'^'^ and H'^'P' denote the planar and non-planar parts, respectively. In comparison 
to the case of a real scalar field, the numerical prefactors of planar and non-planar parts 
have changed. 

The integral (87) has been evaluated in reference |32| by using Schwinger's exponential 

<k 

P 



parametrization, the decomposition cos{kp) = 2 S e^'''^^ and 

r]=±l 



+ m2 + I2 



M4 



2M2 



2 ^ fc2 I 



+ CM2 



(88) 



®In order to simplify computations, we consider the matrix (S*"^) to be block-diagonal so that — 9^k^, 
see Eqns. ©jS. 
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where 



(which may be real or purely imaginary depending on the value of 



the parameter a). For p 1, the non-planar part behaves like 



n-pi(p) 



4(47r) 



p2 



+ m In ( p 



In 



M2 



+ o{i), 

(89) 



and thereby involves a quadratic IR divergence (and a subleading logarithmic IR diver- 
gence). For a — )■ (i.e. — )• ^) this result reduces to the one for the naive model, 
i.e. 



lim n"-p'(p) 

a— !>0 



-A, 



4(47r)^ 



p2 



+ In (p^m^) 



+ 0(1). 



(90) 



The integral defining the planar part does not contain a phase factor and is therefore 
UV divergent. It can be regularized by introducing a cutoff A and subsequently taking the 
limit p2 _). 0, as explained in reference |32j. The final result can be expanded for large 
values of A, yielding 



rcgul. 



(A) 



-(A1 + A2) 
4(47r)2 



4A2 + m2 In ( ^ 

M2 + 



4M2 



In 



-hM2 
-M2 



+ 0(1). 



(91) 



The one-loop renormalization can be discussed following the case of a real scalar field, 
see reference or section 4.2, Accordingly, the action (85) is expected to be stable with 



respect to radiative corrections to all orders. 

We now turn to the symmetries of the theory. The invariance of the tree level action 
(85) under rigid C/(l)-transformations is expressed by the classical Ward identity 

6 6 



Wr(°)[y?,(^] =0, with W 



6(p 6(p 



(92) 



The tools of algebraic renormalization allow us to investigate whether this symmetry holds 
to all orders of perturbation theory. For this investigation we use the following notation 
for the loop expansion of the vertex functional: 



n-l 



(n-l) 



(m) 



(93) 



m=0 



Using the Ward operator W and the charge conjugation C : ip (p, we now try to show 
recursively that the radiative corrections will not destroy the rigid U{1) symmetry of the 
vertex functional at higher order. At order (n — 1) in ^, we have 



(94) 



where we assumed that we can ensure the symmetry with appropriate counterterms at 
the (n — 1) loop order. In the last equality of Eqn. (94) we have assumed the validity of 
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the QAP: the symmetry breaking term A then represents a Poincare invariant, integrated 
insertion with a mass-dimension less or equal four. In contrast to the commutative case 
(where A must be local) , we presently allow for non-local insertions of the type discussed 
in Section [5} We note that the expression A is odd with respect to charge conjugation (i.e. 
CA = —A) as a consequence of the fact that the Ward operator W is also odd under the 
operation C. 

The classical field functionals A obeying all of the constraints we just listed can be 
expanded with respect to the basis given in Eqn. (79). Obviously, the five basis elements 
can be written as the W-variations of the following five independent field functionals A: 



1 1 d'x {^^'^ + 



d X LD — LD + Lp — ip 
□ □ 



Thus, one has 



WA,: = A,: 



(95) 



(96) 



which implies that we can rewrite (94) as 



with A 



1=1 



(97) 



Let ^(^.i) denote the action which includes all counterterms up to order (n — 1) and which 
leads to the vertex functional T(^n-i)- If we now replace the action 8(^^-1) by tlie new action 



in) 



(n-1) 



then we get the new vertex functional 



in) 



(n-1) 



h'^A + 0{K 



n+l\ 



SO that (94) yields 



wr(„) = + o{K' 



n+l\ 



(98) 
(99) 
(100) 



Thus, we have verified recursively that the Ward identity holds to all orders and that there 
is no anomalous breaking of symmetry. 

We stress that the previous procedure can only be expected to work for models for 
which one has taken care of the UV/IR mixing problem. Thus, for the considered example 
it is crucial to include the 1/p^-term in the action which implements the IR damping in 
the propagator. 



7 Conclusions and final remarks 

In this paper we have shown the validity of the QAP at the lowest orders of perturbation 
theory for the case of a non-commutative scalar field theory. Although we expect the 
QAP to hold in general for renormalizable models, a complete proof still remains to be 
elaborated. 
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In the sequel, we have discussed the difficulties that one encounters when applying the 
usual BPHZ approach to non-commutative field theories. In the commutative case, the 
usual BPHZ subtraction scheme together with Zimmermann's forest formula [38| itera- 
tively solves all problems raised by multi-loop graphs with overlapping divergences. For 
the non-commutative case, we have put forward a procedure to remedy these problems 
and we have explicitly shown that our strategy works at the one-loop level. We are cur- 
rently investigating the so-called sunset graph which is a typical two-loop graph with an 
overlapping divergence [S^. Indeed overlapping divergences represent a delicate issue in 
renormalization which have to be dealt with in any scheme. 

Finally, we have analyzed the nature of non-localities which can appear in typical non- 
commutative field theories and we applied the method of algebraic renormalization while 
taking into account these non-localities. We hope that this program can be further com- 
pleted so as to be applied to non-commutative gauge field theories whose renormalizability 
represents an outstanding problem. 
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A Evaluation of expression (50) 



To evaluate expression ( 50 ) , we write 



:ioi) 



ri=±l 



and we use Schwinger's parametrization: 

2kp + 



J{p) = - (fk 



{{k + pY + m'^){k'^ + m^)^ 



cos(A:^) = - ^ Jr,{p) 



(102) 



ri=±l 



with 



oo oo 



J^(p) = j (fk {2kp + p'^) jdaa j dj3 exp [-a [k"^ + vn?) - I3[{k+ pf + m^] + ir]kp] . 



(103) 

In order to carry out the A;-integral in Jrjip), we complete the square in the exponent: 



oo oo 



Jrj{p) = J daa j dp j d k' {2kp + p j exp 





-(« + /3)fc" + ^^^^-/3/-(a + /3)m2 

(104) 
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where 



k' = k + 



Q + /3 



(105) 



After expressing the factor {2kp + p^) in terms of k' , using pp = 0, and performing the 
Gaussian integration over k\ we get 



oo oo 



Jrf{p) = —TT^p'^ J da J dp 





a 



2/3 

a + /3 ^J{a + PY 



exp 



a + P 



Pp^ - (a + P)m^ 



The change of variables (a,/3) — )• (A,^) defined by 

a = (l-^)A, P = i\, with A G [0,oo[, ^£[0,1], 
hence a + /3 = A , dadp = XdXd^ , 

leads to 



(106) 



(107) 



Up) = -^V /^^ (1 - 26(e - 1) e-'^^Pf y^A exp - A [m^ - (^ - 1)^^' 



(108) 



where we used r/^ = 1, and by virtue of = we have e "'^^^ = — 1. Carrying 



77=±1 



out the integration over A we end up with an integral over a modified Bessel function Ki{z) 
of second kind for expression (102): 



Jip) = 7:V I 



1 (1 - 2C)(e - 1) K, ( v(m2 - (e - i)ep2)p2 



(109) 
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